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ABSTRACT 

Strong size and internal density evolution of early-type galaxies between z ~ 2 and the present has 
been reported by several authors. Here we analyze samples of nearby and distant [z ~ 1) galaxies 
with dynamically measured masses in order to confirm the previous, model-dependent results and 
constrain the uncertainties that may play a role. Velocity dispersion (cr) measurements are taken 
from the literature for 50 morphologically selected 0.8 < z < 1.2 field and cluster early-type galaxies 
with typical masses Md yn = 2 x 10 11 M@. Sizes (R c s) are determined with Advanced Camera for 
Surveys imaging. We compare the distant sample with a large sample of nearby (0.04 < z < 0.08) 
early-type galaxies extracted from the Sloan Digital Sky Survey for which we determine sizes, masses, 
and densities in a consistent manner, using simulations to quantify systematic differences between 
the size measurements of nearby and distant galaxies. We find a highly significant difference between 
the a-R c ff distributions of the nearby and distant samples, regardless of sample selection effects. The 
implied evolution in R e g at fixed mass between z = 1 and the present is a factor of 1.97 ± 0.15. This 
is in qualitative agreement with scmianalytic models; however, the observed evolution is much faster 
than the predicted evolution. Our results reinforce and are quantitatively consistent with previous, 
photometric studies that found size evolution of up to a factor of 5 since z ~ 2. A combination of 
structural evolution of individual galaxies through the accretion of companions and the continuous 
formation of early-type galaxies through increasingly gas-poor mergers is one plausible explanation of 
the observations. 

Subject headings: galaxies: clusters: general — galaxies: elliptical and lenticular, cD — galaxies: 
evolution — galaxies: formation — galaxies: fundamental parameters — galaxies: 
general — galaxies: photometry 



1. INTRODUCTION 

Hierarchical galaxy formation models embedded in a 
ACDM cosmology predict strong size evolution for mas- 
sive galaxies. A higher gas fraction in high-redshift galax- 
ies leads to more dissipation and hence compact galax - 
ies fe.g.. lRobertson et al.ll200llKhochfar fc SiUdl2006bD . 
and subsequent evolution such as dry merging or accre- 
tion of smaller systems can increase the size of a galaxy 
(e.g. lLoeb fc Peebledl2003l iNaab et al J 12001 . The mod- 
els predict the strongest sampl e-averaged size evolution 
for the most massive galaxies (|Khochfar fc Silk) l2006aT> 
because of large differences in the gas fraction at different 
redshifts and because the assembly of massive galaxies 
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continues until very late epoch s in a hierarchical frame- 
work (e . g ., EeL^ciairin |20Q3) . 

Evidence for significant size evolution between z ~ 
2 and the present has be en building up quic k ly over 
the past few ye ars (e.g. iTruiillo et alt 12004 l2006bt 
iFranx et al.ll2008f) . In particular, galaxies with low star 
formation rates and high stellar masses (> 1 11 ap- 
pear t o be extremely comp a ct from z ~ 1.5 (|Daddi et al.l 
120051: ITruiillo et all l2006al l2007t iLonghetti et all 120071 : 
iCimatti et all 120081: iRettura et all 120081) to z- ~ 2.5 



(IZirm et al.l l2007t iToft et alj l2007t Ivan Dokkum et all 
120081: iBuitrago et al.ll2008f) . Given the similarity between 
many of their observed properties there is likely to be an 
evolutionary connection between these distant compact 
galaxies and the present-day early-type galaxies despite 
the measured large differe nce of 2 orders of m agnitude in 
surface mass density (e.g.,[Zirm eFaT]l200l . 

The measurements of sizes and densities of high- 
redshift galaxies are hampered by many systematic un- 
certainties (e.g., morphological if-corrections, surface 
brightness dimming, errors in photometric redshifts, and 
mass measurements). Most of these errors, however, are 
unlikely to fully account for the observed strong size evo- 
lution. The uncertainty in the mass estimates may be the 
exception. For the work in the literature, these mass es- 
timates are always based on the photometric properties 
of the galaxies. For a reasonable set of assumptions the 
photometric stellar mass estimates are not uncertain by 
more than a factor of 2 or 3 and would not change the 
inferred evolution significantly. However, since we infer 
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that z ~ 2 galaxies must have physical central densities 
that are 3 orders of magnitude hig her than those of local 
galaxies (|van Dokkum et alll2008[ ). further verification of 
those apparently reasonable assumptions is warranted. 
For example, a stellar in itial mass function (IMF) that is 
radically different (e.g.. iLarsonl 120051 : iFardal et all 120071 : 
Ivan Dokkumll2q08t lDave!l2q08l) from a Salpeter-likg IMF 
dSarpeterfll955l IScaldll986t lKroupal l200ll; IChabrienl2003l; 
IHoversten fe Glazebrook! I2008D could reduce the stellar 
mass estimates by an order of magnitude, producing per- 
fectly normal galaxies by today's standards. 

The spectacular nature of these compact galaxies 
at z ~ 2 could be confirmed by direct, kinemati- 
cal mass measurements. However, the quiescent na- 
ture of th ese objects and th eir consequent lack of emis- 
sion lines (|Kriek et al. 2006) require absorption-line mea- 
surements of their stellar velocity dispersions, which 
should be as high as 400 - 500 km s _1 (jToft et all 
120071: Ivan Dokkum et aTll2008l ). Unfortunately, with the 
currently available instrumentation this is not feasible. 
These z ~ 2 galaxies are proh ibitively faint at obs erved 
optical wavelengths (see, e.g.. ICimatti et all 12008). and 
near-infrared spectroscopy is still maturing as a tech- 
nique. Continuum detections i n the observed N IR have 
only recently become possible (jKriek et al.ll2006T ) for the 
brightest sources, and no detection of absorption lines 
has been made. 



At lower redshifts [z ~ 1) absorption- line spectroscopy 
has for years been a powerful t ool to study the evolution 
of distant early-type gala x ies | van Dokkum et "ail 119981: 
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process (see, e.g.. ITruiillo et all l2006bl ): therefore, 
intermediate changes in sizes and densities should be 
observable at these redshifts. 



In this paper we compile a sample of galaxies at red- 
shifts 0.8 < z < 1.2 with measured absorption line veloc- 
ity dispersions from the literature and that are visually 
classified as early-type galaxies with the aid of Hubble 
Space Telescope (HST) i maging from the Advanced Cam- 
era for Surveys (ACS; iFord et a l. 1998). We measure 
the galaxies' sizes from these ACS data. We then com- 
pare this distant sample with nearby early-type galax- 
ies extracted from the Sloan Digital Sky Survey (SDSS; 
I York et al.ll200l . This comparison, with careful con- 
trol of systematic uncertainties, allows us to verify that 
distant early-type galaxies are indeed significantly more 
compact than their local counterparts. 

The advantage of this approach is that the size and 
density measurements are independent of the photomet- 
ric properties of the galaxies apart from the surface 
brightness profile. The absence of luminosity and other 
photometric properties from our analysis assures us that 
our study does not suffer from the strong possible biases 
in previous photometric work. Moreover, deep, high- 
resolution ACS imaging allows us to determine sizes of 
z ~ 1 galaxies to a precision comparable to that nearby 
galaxies and in a consistent manner. Most previous stud- 
ies verify for biases in the size determinations within 
their distant samples (e.g.. ITruiillo et all l200i |2006bt 



ICimatti et alll2008l ) but do not extend this analysis to 
verify the consistency with size measurements of nearby 
galaxies. 

In § [2] we describe the sample of nearby early-type 
galaxies and derive the dynamical mass-size relation. In 
§[3] we construct the sample of z ~ 1 early- type galaxies, 
determining their masses and sizes in a manner that is 
consistent with the nearby sample. In § [4] we quantify 
systematic effects in our size measurements through sim- 
ulations. In §[5] we derive the evolution in the dynamical 
mass-size relation. In § [6] we compare our results with 
previous measurements based on photometric mass esti- 
mates and semianalytic model predictions, and we dis- 
cuss size evolution in the broader context of the evolving 
early-type galaxy population. Finally, in §[7] we summa- 
rize our results and conclusions. We adopt the following 
cosmological parameters: (Qm, ^a, h) — (0.3, 0.7, 0.7). 

2. NEARBY EARLY-TYPE GALAXIES 

2.1. Velocity Dispersions and Sizes 

We have extracted a large sample of early-type galax- 
ies at redshifts 0.04 < z < 0.08 from t he SDSS database 
(DR6; lAdelman-McCarthv et alll2008l) based on the cri- 
teria as outlined bv lGraves et al.l (|2007f ). 5 Briefly, galax- 
ies on the red sequence and either without emission lines 
or with high [Oil] to Ha ratios are included in the sam- 
ple. These criteria effectively exclude star-forming galax- 
ies, but incl ude genuine early -type galaxies with nuclear 
activity (see lYan et alll2006D . 

The dispersion as measured within the spectroscopic 
aperture (cr a p) is corrected to match the average disper- 
sion within the effectiv e radius R e g (measure d as de- 
scribed below) following I Jorgensen et al.l (|1995f ): 



where R ap (z) is the radius of the SDSS spectroscopic 
fiber (1.5") in kpc at t he distance of the ga laxy. We 
use the correction from Jorgcn sen et all (| 19951 ) for con- 
sistency with previousl y published results. We note that 
Cappellar i et al.l (|2006f ) used better data to improve the 
aperture correction, but the resulting difference in <j e g is 
only a few percent. 

We use GA1FIT (jPeng et all 120021 ) to determine ef- 
fective radii from the SDSS g-band imaging assuming 
an i? 1 / 4 profile, leaving the effective radius, the inte- 
grated magnitude, the position angle, the axial ratio, 
and the position of the center as free parameters. The 
point-spread function (PSF), which is used to decon- 
volve the image, is constructed for each galaxy separately 
by co-adding the stars in the frames after drizzling the 
cuto uts to a com mon center. A more general R}l n pro- 
file (|Sersid Il968l) may provide a more realistic descrip- 
tion of the surface brightness distribution of individual 
early- type galaxies, especially in the presence of a signif- 
icant disk. However, n = 4 provides a good description 
of the average profile of early-type galaxies both nearby 
(|de Vaucouleurslll948f ) and at z ~ 1 (see § EOj) . More- 
over, introducing n as an additional free parameter re- 
sults in unnecessarily large, redshift-dependent system- 
atic uncertainties in the size measurements (see § |4|) . 

5 The IDs, positions, redshifts, and dispersions for this sample 
were kindly provided by G. Graves 
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The size parameter that we use in this paper is the 
circularized effective radius y/ab = a^q, where a is the 
effective radius along the major axis (the output param- 
eter of GALFIT, b is the effective radius along the minor 
axis, and q the axis ratio (as calculated by GALFIT); 
\fab is a good approximation for optically thin luminos- 
ity distributions such as the generally dust-poor early- 
type galaxies in our samples. The systematic and ran- 
dom errors of our size determinations are inferred from 
extensive simulations described in § 21 

The SDSS spectroscopic catalog suffers from several 
biases that may mitigate size evolution measurements. 
First, compact sources are not targeted for spectroscopy 
as they may be mistaken for stars or because their central 
surface brightnesses, i.e., their fiber magnitudes, are too 
bright. Second, almost all galaxies that in the literature 
(see t he HyperLEDA database compiled bv lPaturel et al.l 
2003;) have been claimed to have high, > 300 km s -1 ve- 
locity dispersi ons have disper sions of < 300 km s -1 in 
the SDSS fsee lBernardil [20071 Appendix A). The source 
of this discrepancy is unknown. While Bernardi convinc- 
ingly argues that the SDSS dispersions are more reli- 
able, there are a number of galaxies with large, mutually 
consistent dispersion measurements from multiple, inde- 
pendent observers, and for which the SDSS dispersion is 
significantly smaller. These potential caveats may cause 
our size evolution measurements to be biased. We refer 
to these issues when we present our results in § [5j 

2.2. The Mass-Radius Relation and the Mass-Density 

Relation 

From er ff and R c g we derive the total dynamical 
mass and the corresponding average surface mass den- 
sity within R c g : 
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ies (|Cappellari et al.ll2006f) . Following IShen et all (2003) 
we adopt the following characterization of the Md yn --R e ff 
relation: 

With a least-squares linear fit to all galaxies with mass 
Af dyn > 3 x 10 10 M Q we find that the slope is b = 0.56 
and the zero point normalized to a characteristic mass 
M c = 2x 10 11 M Q is R c = 4.80 kpc. We find statistically 
the same relation if we perform a linear fit to the values 
of the median R g in 0.1 dex wide mass bins in the range 
10.5 < log(Mdyn) < 12.1. The scatter around the best-fit 
relation is 0.14 dex. 

Using stellar masses, M„, IShen et, all (l2003h find the 
same slope b = 0.56 for the M*-R c g relation. Their 
zeropoint, however, is larger (R c = 6.14 kpc). This is 
likely due to th e diffe rence betwe en Md vn and M „ as 
ICappellari etafl (|2006h show for a iKroupal (|2001h IMF 
(which is also used by Shen et al. 2003) that Md yn ~ 
1.4M*. This translates into a difference in R 9 of ~20%, 
close t o the observed difference. Furthermore. IShen et al.l 
(20031) analyze SDSS r-band imaging whereas we use g- 
band imaging. 



3. DISTANT EARLY- TYPE GALAXIES 
3.1. Velocity Dispersions and Sizes 

Several authors have published velocity dispersion 
measurements of early- type galaxies at z ~ 1. We com- 
pile the data from three different data sets for which 
the selection criteria are well understood so that sys- 
tematic effects introduced through selection effects can 
be properly modeled. Our compiled sample contains 
galaxies in the redshift range 0.8 < z < 1.2 in the 
Chandra D eep Field-South (CDF-S, Ivan der Wel"et~all 
1 2004 12005h and the Hubble Deep Field-North (HDF-N, 
ITreu et all l2005allbr ) . In addition, we include galaxies 
in the mas sive, X-ray selecte d cluster MS 1054-0321 at 
z = 0.831 ((Wuvts et all 12004( 1. The seven cluster galax- 
ies at z > 1 for which dispersions have been meas ured 
(jvan Dokkum fc Stanford! l2003t iHolden et all 12005( 1 are 
not included because of the paucity of this sample, which 
prevents us to accurately model selection effects. The fi- 
nal sample only contains galaxies with S/N > 10 A 
since dispersions derived from spectra with lesser quality 
can suffer from large (> 10%) systematic uncertainties. 
The same aperture corrections are included as for the 
nearby galaxies (Eq. []}, with in this case R ap (z) the ra- 
dius of a circle with area 1" x 1.25" (the width of the 
slits and the typical length of the extracted region) in 
kiloparsec at the redshift of the galaxy. The data are 
given in Table [TJ 

For all galaxies ACS imaging is available. In order to 
produce an internally consistent data set, we re-measure 
the sizes of all galaxies. GOODS 6 provides deep, pub- 
licl y available ACS imagin g of the CDF-S and the HDF- 
N (jGiavalisco et all 12004(1 in four filters. We use the 
F850LP (z 850 -band) images in order to match the rest- 
frame wavelength at which the sizes of the nearby com- 
parison sample are measured (the SDSS g-band; see §[2]). 
For the MS 1054-0321 cluster ACS imaging has been 
taken as part of the guara nteed time observation pro- 
gram (jBlakeslee et al.1 12006( 1. We use the F775W (i 775 - 
band) imaging as the available z-band imaging is of lesser 
quality. At this redshift rest frame g falls in between 
1775 an d -2850 such that the morphological if-correction 
is not a problem; the z ~ 0.8 galaxies in the sample of 
ITreu et al.l (|2005bl l are only 3% ± 4% smaller in the i 775 
band than in the zsso band. 
With GALFIT we determine effective radii in the same 
for the nearby galaxy sample (see jj |2~Tj) . 
constructed with Tiny Tim (|Kristl )l995). 
even thou gh using stars resu l ts in negligible di f ference s 
(see, e.g.. Ivan der Wei et al.l 120051 : ITreu et all l2005bt l. 
Errors are discussed in § 21 and the data are given in 
Table [Q 

We use Eqs. [2] and [3] to compute masses and sur- 
face densities. Again, we adopt = 5, which 
has been shown to hold for distant nonrotat i ng el- 
liptical galaxies (Ivan der Marel fc van Dokkuml l2007t 
Ivan der Wei fc van der Marel 12008( 1. For rotating early- 
type galaxies the situation appears to be more complex 
(|van der Wei fc van der Marelll2008( ) in the sense that j3 
is possibly ~20% larger than 5. We comment on the im- 
pact of this possible complication on our size evolution 
measurement in S 15.21 A low-mass cutoff of 3 x 10 10 Mr?, is 



manner as 
The PSF is 



See http://www.stsci.edu/science/goods/ 
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Fig. 1. — Stacked zssrj-band image of 29 z ~ 1 elliptical galaxies 
in the CDF-S and the HDF-N and the residuals after subtracting 
J? 1 /" model profiles, with n = 3, 4, 5, 6, 7. Each panel is 7.68" on a 
side, which corresponds to 62 kpc at z = 1. The circle indicates the 
effective radius as measured with the i? 1 / 4 model profile (0.31", or 
2.5 kpc at z = 1). The model with n = 4 provides the best fit. The 
model with n = 3 produces positive residuals at large radii; models 
with high n produce negative residuals. This justifies our choice 
to adopt the i? 1 / 4 law to model the surface brightness profiles of 
both local and distant early-type galaxies. 

applied since below this limit no useful samples are avail- 
able due to severe incom pleteness of the surveys (see, e.g., 
Ivan der Welli~all[2005h . 

3.2. The Average Surface Brightness Profile 

In determining the sizes of the nearby and distant 
galaxies in the previous sections we assumed that an 
i? 1 / 4 profile provides an accurate description of early- 
type galaxies. We know this to be true for nearby galax- 
ies, but not for z *~ 1 early-type galaxies. If a more 
general R}l n profile is adopted, measured values tend 
to cluster around n = 4 (see, e.g.. iBlakeslee et al.ll2006t 
iRettura et al-l feOOG) . However, there is a possibility that 
the true values of n are different; at large radii the 
"wings" of the profile become quickly overwhelmed by 
background noise, even in the deepest HST imaging. Be- 
cause n and the measured R c g are correlated, assuming 
n = 4 for all redshifts introduces systematic errors in 
case n evolves with redshift. 

In order to examine the profiles of the z ~ 1 galax- 
ies at large radii, we median-stack the zs50 _ band images 
of all elliptical galaxies (SOs are excluded) without bright 
neighbors in our CDF-S and HDF-N samples (see Fig. Q]) . 
The images of the individual galaxies are drizzled onto 
a common central position. Due to imperfections in this 
procedure, the stacked PSF may not be an accurate de- 
scription of the PSF of the stacked image. However, this 
does not play a role since the deviations we are interested 
in have scales that are an order of magnitude larger than 
the PSF. 

With GALFIT we subtract R 1/n profiles with integer 
values n = 3 — 7 (see Fig.[lJ. The negative residuals out- 
side i? c ff for models with large n and the positive resid- 
uals for models with small n indicate that these limiting 
cases provide poor fits of the outer regions of elliptical 
galaxies at z ~ 1. The i? 1 / 4 and i? 1 / 5 profiles provide 
the best description of their average surface brightness 
distributions. This visual impression is confirmed by the 
X 2 -values of the respective fits: x 2 = 0.5 for both n = 4 
and for n = 5, whereas \ 2 > 0.7 for other values of n. In- 



terestingly, n does not evolve significantly with redshift, 
and we conclude that it is safe to assume that choosing 
n = 4 for both nearby and distant early-type galaxies 
does not introduce significant systematic errors. 

4. SIMULATIONS OF SIZE MEASUREMENTS 

To test the robustness of our size determinations of lo- 
cal and distant early- type galaxies in § § [2] and [3] we sim- 
ulate size measurements by using SDSS q-ban d imaging 
of 45 early-type galaxies in the Virgo Cluster (|Mei et al.l 
l2007h . The pixels of the mosaics of the Virgo Cluster 
galaxies are re-binned to account for the different pixel 
scales of the various instruments, and different cosmo- 
logical distances of the galaxies at higher redshifts. The 
redshift range is z — 0.04 — 0.08 for the nearby sample 
and 0.8 — 1.2 for the distant sample. The physical sizes 
of the simulated galaxies are thus conserved. For each 
redshift [z = 0.04,0.06,0.08,0.8,1.2) we run -200 sim- 
ulations with different values for the flux density of the 
simulated galaxies, which are chosen such that the sim- 
ulated galaxies have the same range in apparent mag- 
nitude as the observed galaxies in our samples. After 
convolution with the appropriate PSF and the addition 
of Poisson noise the seed galaxies are inserted into empty 
parts of real images. Their sizes are measured with GAL- 
FIT in the same manner, by fitting a i? 1 / 4 law, as for the 
real galaxies. Because we are mainly interested in the 
systematic differences in the size determinations within 
and between our nearby and distant samples, we assume 
the size determinations based on the z — 0.04 simulated 
SDSS images of the Virgo Cluster galaxies as the base- 
line against which we compare the other simulated size 
measurements . 

The results of the simulations are shown in Fig. [5J 
Within the nearby sample we find systematic, redshift- 
dependent differences, of up to ~10%. Random errors, 
derived from the scatter in the sizes inferred from the 
simulated images, are typically less than 5%. System- 
atic differences between the nearby and distant samples 
can be as large as 20% for small galaxies, where at high 
redshift the sizes are overestimated. Random errors are 
typically 10-15%. We find no systematic trends with 
magnitude. The reason for this is that all galaxies are 
relatively bright compared to the depth of the data sets, 
such that the limiting factor in the size measurements is 
spatial resolution. 

Adopting a R}l n law with n as a free parameter may 
improve the quality of the fits to individual galaxies. 
However, our simulations reveal that the random errors 
increase to ^20-25%, without much change in the sys- 
tematic errors. Together with the analysis of stacked 
images (§ 13. 2| . this test justifies our choice to use the 
i? 1 / 4 law to describe the surface brightness profiles of all 
galaxies in both the nearby and the distant samples. 

The sizes we derive in the §§[5] and [31 and the derived 
quantities Md yn (Eq. [2]) and £ e ff (Eq. [3]) , are corrected to 
account for systematic measurement errors. Those cor- 
rections depend on redshift and are interpolated between 
the values listed in Fig. [51 For simplicity the dependence 
on size is not taken into account, such that the remaining 
systematic uncertainty is ~5%. 

5. STRUCTURAL EVOLUTION OF EARLY-TYPE 
GALAXIES 
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Fig. 2. — Random and systematic errors in the size determi- 
nations of early- type galaxies with SDSS imaging at (a) z = 0.06 
and (6) z = 0.08 and with HST imaging at (c) z = 0.8 and (d) 
z = 1.2, all with respect to the size measurements at = 0.04, which 
are used as the benchmark in our analysis. These are the results of 
simulations with 45 early-type galaxies in the Virgo Cluster. The 
systematic offsets are indicated by the dashed lines. The scatter in 
the offsets, considered to be the random error in the size determi- 
nations, are also listed. 

5.1. Evolution of the a-Radius Distribution 

In Figure [3] we compare the <T e ff--Reri distributions 
of the nearby and distant early-type galaxy samples. 
This unusual projection of the fundamental plane (FP; 
iDressler et aLlll987tlDiorgovski fc Davislll987t ) has a very 
large scatter. However, the advantage is that changes 
with redshift are independent of luminosity evolution. 
Despite the large scatter, it is clear that the distant 
galaxies are offset from the nearby galaxies. Galaxies 
with the properties of typical galaxies in the distant 
sample (a c s ~ 250 kms -1 ; R c g ~ 3 kpc) are rare 
in the local universe. In the nearby sample, galaxies 
with (T c ff ~ 250 kms -1 have much larger sizes, and 
galaxies with sizes R e g ~ 3 kpc have dispersions of 
Ccff ~ 150 kms -1 . These numbers are only intended to 
guide the eye. A quantitative analysis of the differences 
between nearby and distant galaxies is presented below. 

The distant sample is not directly comparable with 
the nearby sample in its entirety (left), as the nearby 
sample reaches to much lower masses. In order to as- 
sess the question whether the true, underlying a c {f-R c g 
distribution of distant galaxies is different from the a D g- 
i? e ff distribution of nearby galaxies, we need to remove 
the galaxies in the nearby sample that would not be in- 
cluded at z ~ 1 due to sample selection effects. The 
sub-sample of nearby galaxies that is observable at z = 1 
is shown in the right panel of Fig. [3] The two criteria 
that the galaxies in the observable sub-sample satisfy are 
L > L min and i£ e ff < ^efr.max- L mm ~ 10 10 L Qi b is the 
luminosity limit for the field z ~ 1 surveys (see, e.g., 
Ivan der Wei et aT] 120051 ) after correcting for luminosity 
evolution between z = 1 and the present (0.555 dex; 



Ivan D okkum fc van der Marej|2007ri. Fo r the M S 1054- 
0321 cluster sample from |Wivts_et_alJ (|2004f ) this is 
1.8 x 10 10 i Q! B- The second criterion i? ff < i? e ff,max 
takes into account that high signal-to-noise ratio (S/N) 
spectra are harder to obtain for low surface brightness 
galaxies than for high surface brightness galaxies with 
the same luminosity; i.e., the distant sample is biased 
in favor of small ga l axies. Th e S/N of t he spe ctra 
of Ivan der Wei et all (|2005l ) and iTreu et "all (|2005bh do 
not precisely scale linearly with luminosity L = IR 2 S , 
where I is the surface brightness, but as S/N oc IR\§- 
This implies that, at fixed luminosity L, S/N oc R~ 0A . 
Since a dispersion measurement requires a minimum S/N 
(~ 12A" 1 ), a galaxy with luminosity L has a maximum 
radius i? e ff,max oc L 2 5 for which its dispersion can be 
determined. We use the luminosity limits of the surveys 
discussed above to normalize the dependence between lu- 
minosity and maximum size; we simply assume that for 
the smallest galaxies (i? c ff = 1 kpc) the luminosity limit 
coincides with the size limit such that we have 

■Rcff.max (kpc) = ( ) . (5) 



One would expect that for galaxies smaller than 3 kpc 
the signal-to-noise ratio of the spectra would not depend 
on size any longer since seeing generally dominates the 
apparent sizes of such small galaxies at z ~ 1. Because 
of the variety of telescopes, weather conditions, and data 
reduction techniques, this, however, is washed out and 
not apparent in the data. We note that the introduction 
of, effectively, a rudimentary surface brightness criterion 
is a step forward in modeling the selection effects with 
respect to earlier attempts that only take total luminosity 
into account. 

The difference between the ir ff-i? ff distributions at 
low and high redshift is highly significant, even after tak- 
ing selection effects into account (Fig.[3J right). The two- 
dimensional Kolgomorov-Smirnov statistic has a high 
value (D = 3.71), which implies that it is extremely un- 
likely that the nearby and distant samples are drawn 
from the same distribution. By repeatedly drawing sam- 
ples from the nearby sample with the same size as the 
distant sample we confirm this: less than 0.001% of the 
simulated samples have D = 3.71 or higher. 

5.2. Evolution of the Mass-Radius Relation 

The structural difference between the nearby and dis- 
tant samples described in the previous section implies 
that the Md yn -i? c ff and Md yn -S e ff relations evolve with 
redshift. In Fig. [4] we show the -Md yn --Rcff relation for the 
distant sample, and compare this with the equivalent re- 
lation for nearby galaxies derived in § 12.21 Clearly, the 
relation shifts to smaller radii from low to high redshift. 

Parametrized as in Eq. [4] we find R c = 2.58 ± 0.17 and 
b = 0.65 ± 0.06 with a scatter of 0.117 ± 0.013 dex (af- 
ter subtracting the observational uncertainties in quadra- 
ture). The errors are estimated with a bootstrap/Monte- 
Carlo simulation in which the data points are randomly 
sampled and varied according to the (correlated) mea- 
surement errors which are assumed to be Gaussian. The 
systematic error in i? ff of 5% (see § [4j) is also taken into 
account. 

The treatment of the selection effects described in the 
previous section shows that the observed size evolution 
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Fig. 3. — The cr e ff-i? of j distributions of the nearby sample of early-type galaxies (gray scale) and the distant sample (data points). 
The red data points are cluster galaxies, and the blue data points are field galaxies. The error bars at the top right indicate the typical 
values of the errors for the distant sample. The solid line indicates log(Md yn ) = 10.5. Dotted lines indicate lines of constant M^ yn (parallel 
to the solid line), spaced by 0.5 dex, and lines of constant surface density £ e ff (perpendicular to the solid line), also spaced by 0.5 dex. 
The left-hand panel shows the entire nearby sample; the right-hand panel only shows those g alax ies in the nearby sample that would be 
included in the distant sample considering the selection effects, that apply to the surveys (see § 15. 11 1. The highly significant offset (> 99.9% 
significance) between the local and distant samples implies significant structural evolution in the early-type galaxy population between 
z ~ 1 and the present. 



seen in Fig. [5] is not an artifact. However, given the na- 
ture of the selection effects, which favor small galaxies 
over large galaxies, the intrinsic amount of size evolu- 
tion and possible evolution in the slope and scatter of 
the -Md yn --Rcff relation must be inferred through careful 
modeling. The goal is to derive the intrinsic Md yn --Reff 
relation at z ~ 1 that reproduces the observed M ( j yn -i? e ff 
distribution after applying the selection criteria. We take 
an iterative approach due to the interdependence of the 
selection criteria and the amount of evolution in zero- 
point, slope, and scatter of the Md yn --Reff relation. In 
the following we de-evolve the properties of the nearby 
sample to constrain the form of the true, underlying 
z ~ 1 Ma yn -R e ff relation. 

The simplest evolutionary scenario is a change in the 
zero point R c (see Eq. 2J. For each object in the 
nearby sample the size is reduced by the same amount 
Alog(-Reff), and those that do not satisfy the selection 
criteria described in the previous section are removed. 
From the remaining sub-sample the "observed" Md yn - 
R c e relation is determined. The different selection cri- 
teria and sample sizes for field and cluster galaxies are 
taken into account in this process, which is repeated for 
many different values of Alog(i? c ff)- We find that an 
intrinsic value of R c = 2.64 ± 0.18 reproduces the ob- 
served value of R c = 2.58 ± 0.17. Hence, it appears that 
selection effects do not strongly affect the inferred size 
evolution. 

However, the scatter of the assumed intrinsic distri- 
bution (0.14 dex) is higher than the observed scatter 
(0.117 ± 0.013 dex). This cannot be explained by se- 
lection effects in the simple scenario described above. It 
is therefore required that the scatter, as well as the zero 
point, is also treated as an evolving parameter. This is 



implemented in our analysis by reducing or increasing 
the offset of each galaxy in the nearby sample from the 
best-fit -Md yn --Rcff relation by a given fraction. Doing 
so, we find that the best-fitting zero point R c is not dif- 
ferent from the earlier estimate based on a nonevolving 
scatter. We also find that the evidence for evolution in 
the scatter is weak (~ 1.5 a). This exercise mainly serves 
to show that our size-evolution result in not sensitive to 
the amount of evolution in the scatter allowed by the 
observations. 

A similar verification must be carried out for evolution 
in the slope of the Md yn --R ff relation. Allowing only 
the scatter and the size to evolve, as described above, 
the inferred slope of the "observed" Md yn --R e ff relation is 
0.59, marginally consistent with the true observed slope 
of b = 0.65 ± 0.06. If we treat the slope as an addi- 
tional, third free parameter we confirm that evolution in 
the slope, as constrained by our measurements, does not 
affect our size-evolution measurement. An intrinsic slope 
of b = 0.61 provides a better fit than the original slope of 
the Md yn --R e ff relation of the nearby sample (b = 0.56), 
but the difference is marginal (~ 1 a). 

We conclude that, despite (weak) evidence for evolu- 
tion in the slope and the scatter of the Md yn --R e ff rela- 
tion with redshift, there is no significant improvement in 
modeling the observations by adopting slope and scatter 
as free parameters. Modeling the evolution by a frac- 
tional change in size, regardless of mass and offset from 
the local Md yn --Rcff relation, provides an equally good fit. 
Most importantly, changing the slope and scatter within 
the range allowed by the observations does not affect the 
inferred size evolution. We find that R c = 2.64 ± 0.18 
kpc at z = 0.90, a factor of 1.8 ± 0.1 times smaller than 
at z ~ 0.06. 
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Fig. 4. — Mass-size relation for the nearby sample (solid line) 
and at z ~ 1 (dashed line); the symbols are the same as in Fig. [3] 
For the d erivation of the Md yn -R c ft relation for the nearby sample 
see § 12,21 for the d erivation of the M^ yn -R a g relation for the distant 
sample see § 15.21 The smaller, inset panel shows the distribution 
of the two samples around the M^ yn -R c ff relation of the nearby 
sample (the solid line in the large panel). The distant galaxies are 
1.8 ± 0.1 times smaller than the nearby galaxies. It appears that 
the most massive galaxies do not show as large an offset. This 
indicates that size evolution may be slower for the highest mass 
galaxies than for low-mass galaxies, but it has to be kept in mind 
that these very massive galaxies are brightest cluster galaxies and 
may therefore have developed differently from other galaxies. 

The weak evidence for a change in slope of the M^ yn - 
i? e ff relation may also be interpreted as a difference be- 
tween field and cluster galaxies, as the more massive 
galaxies in our sample tend to be cluster galaxies. As- 
suming that slope and scatter remain constant but that 
the zero point of the Md yn --Roff relation evolves differ- 
ently for field and cluster galaxies, we find that R c — 2.49 
kpc for field galaxies and R c = 3.06 kpc for cluster galax- 
ies. The 1 a error on this difference of 0.57 kpc is 0.32 
kpc. The true error may be larger since in this estimate 
it is assumed that scatter and slope behave the same in 
the different environments and that there are no relative 
systematic errors in the size determinations of field and 
cluster galaxies. The evidence for a difference between 
the size evolution of field a nd cluster early-type galaxies 
is therefore weak (see also iRettura et al.l I2008T ) . How- 
ever, we have to keep in mind that so far only a very 
small number of clusters is considered. Future studies 
will need to extend the existing analyses to a larger num- 
ber of clusters to verify the general validity of the results. 

So far, we have assumed that the masses of the galax- 
ies do not change. Our justification is that the scatter 
hardly depends on mass; the effect of a changing mass 
function on modeling selection effects is expected to be 
small. However, physically speaking, it is unnatural to 
propose size evolution without changes in the masses 
of galaxies. Moreover, if the characteristic mass above 
which the number density of galaxies drops off exponen- 
tially evolves with redshift, selection effects will change 
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Fig. 5. — Mass-density relation at z ~ 1. The symbols and lines 
are the same as in Figs. [3] and [4] The % ~ 1 early- type galaxies 
are ~ 4 times more dense than their nearby counterparts. The 
prediction of the semianalytic size-ev olution model for elliptical 
galaxies from Khochfar & Silk ( 2006a) is shown as the dotted line. 
The error bars indicate the predicted size evolution between z = 0.8 
and z = 1.2, the redshift range of our distant sample. Despite 
qualitative agreement, there are significant quantitative differences 
between the predicted and observed evolution. 

as well. The simplest way to implement mass evolution 
is to assume that M cx i? e ff (Eq. [2|) . Including this in our 
modeling procedure shows that the effect on the inferred 
size evolution is less than 5%, and we therefore adopt the 
results with no mass evolution. 

We recall that the nearby sample is biased against 
compact early-type galaxies (§ 12. ip . The potentially un- 
derestimated number of galaxies with dispersions a > 
300 kms -1 is unlikely to drastically affect the size- 
evolution determination for the sample as a whole as the 
average dispersion of the galaxies in the distant sample is 
smaller than that. However, the slope of the local M^ yn - 
i? e ff relation is possibly overestimated, which would lead 
to an underestimate of the slope evolution. More im- 
portant is the problem that small galaxies are missed 
because of their photometric misclassification as stars in 
the SDSS. To fully address this issue a complete analy- 
sis of the SDSS photometric catalog is required, which 
is clearly beyond the scope of this paper. However, we 
can say that it is highly unlikely that the average size of 
nearby early-type galaxies is underestimated by a factor 
of 2 because of this bias. On the other hand, for the 
interpretation of our results and identifying the mecha- 
nisms responsible for size evolution (see § I6.3[) this bias 
could prove to be important. 

In § [3] we noted that the dynamical mass estimate as 
adopted in this paper (Eq.[2j with ,9 = 5) may be too low 
for rotating early-type galaxies. If this is the case, then 
size evolution for these galaxies will be underestimated 
by ~10%. Since this is within the uncertainties of our 
measurements we do not take this further into account. 

Obviously, size evolution at fixed mass translates into 
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density evolution. This is illustrated in Fig. [5] where 
we compare the density distribution of z ~ 1 early-type 
galaxies with the Md yn -S e ff relation for nearby galaxies. 
Because £ e ff does not strongly depend on Afd yn , evolu- 
tionary trends are readily visible; z ~ 1 early-type galax- 
ies are ~4 times more dense than their local counterparts. 
The apparent change in slope can possibly be explained 
by selection effects, completely analogous to our conclu- 
sion that this is the case with the Md yn --R e ff relation. 
Note that compared to the increase in projected density, 
the increase in physical density will be even larger. 

Up until recently, early-type galaxies were thought to 
evolve more or less passively. This appears to be an over- 
simplification and may apply more to their stellar popu- 
lations than to their structural properties. In the follow- 
ing section we discuss possible explanations in the con- 
text of theoretical predictions and the comparison with 
results from studies with different observational strate- 
gies. 

6. DISCUSSION 

6.1. Comparison with Photometric Size- Evolution 
Measurements 

The main goal of this paper is to use dynamical 
measurements to investigate whether early-type galax- 
ies were smaller and denser in the past. Previous work 
has shown that the stellar mass surface density is higher, 
but there are a number of issues with such studies as 
they rely on stellar population models and they ignore 
possible changes in the underlying dark matter profile. 

In Fig. we compare the size-evolution results pre- 
sented in the previous section with size-evolution results 
for early-type galaxies based on photometric mass esti- 
mates. For all the literature samples we take the mean 
rcdshift and the mean stellar mass (normalized to the 
Kroupa IMF), and compute the mean offset from the 
local mass-size relation from lShen et al.l ()2003f ) . We in- 
clude four intermediate redshift cluster galaxy samples 
with photometrically measured masses and sizes from 
WFPC2 or ACS im aging. The data are described by 
iHolden et alJ (|2007h and the sizes are measured as de- 
scribed in this paper (§ [3]). These four clusters are 
CL 1358+62 at z = 0.33, MS 2053-04 at z = 0.59, and 
MS 1054-0321 and RX J0152. 7-1357 , both at z = 0.83. 
Note that we also include MS 1054-0321 in the present 
study with dynamical mass measurements. The agree- 
ment between the independent measurements confirms 
that at least out to z ~ 1 dynamical and photomet- 
ric mass estimates based on optical colors and spectral 
energy distributions agree w ithin the statistica l error s 
as was previously s hown bv Ivan der Wei et alJ (|2006f ). 
iRettura et all (|2006h . and lHolden et all (|2006l ). 

The literature samples have all been selected in dif- 
ferent ways, and so a direct comparison with our work 
may not be straightforward. Not all samples are mor- 
phologically selected; many are selected by their spectral 
or photometric properties. In the local universe there is 
substantial overlap between samples of early-type galax- 
ies that are selected by different criteria; therefore, it is 
a reasonable assumption to suppose this to also be the 
case at high redshift, where different indicators (low star 
formation rates, red colors, smooth visual appearance) 
also reflect a common nature. Recently, several stud- 
ies have shown hints that this is indeed the case (e.g., 



Ivan Dokkum et al.ll2008HKriek et al.ll2008l ). but these is- 
sues need to be further addressed in the future. 

Even with this cautionary proviso, the broad agree- 
ment between the results presented in this paper and 
the photometric results at higher redshifts is striking. 
All studies included in Fig. [6] are consistent with signif- 
icant size evolution of several factors between z ~ 1 — 2 
and the present for galaxies with a given mass. A linear 
fit in log-log space to our two data points at z ~ 0.06 and 
z ~ 1 gives R off (z) tx (1 + ^)-o-98±o.ii With a linear g t 

to the photometric data the inferred rate of evolution is 
i? e ff(z) cx (1 + z )- 1 - 20 ± 012 ; where the error is obtained 
via a bootstrap/Monte Carlo simulation. 

The broad agreement of our measurement of the size 
evolution of early-type galaxies with the photometric 
studies is encouraging and alleviates concerns about se- 
rious systematic effects that potentially could have com- 
promised previous work. Most notably, uncertainties in 
the photometric mass estimates used in all other previous 
work appear to have a limited impact, at least compared 
to factors of >5 which would mimic the strong, observed 
size and density evolution. Uncertainties in photomet- 
ric mass estimates on the level of a factor of ^2 due to 
differ ences among the various s tellar population models 
(e.g.. iBruzual fc Charlotl l200l iMarastonl l2005h remain 
an issue, but to invoke, for example, an unconventional 
stellar IMF as an alternative to radically different struc- 
tural properties of high-redshift early-type galaxies is no 
longer necessary. 

Other systematic uncertainties cannot explain the ob- 
served evolution either. In our size measurements, 
systematic effects have been taken into account (see 
Sees. I3~2l and 14")). We are confident, for example, that 
we would detect low-surface brightness envelopes around 
distant galaxies. Furthermore, we know that only a 
minority of morphologically selected z ~ 1 early-type 
galaxies (~10%) show signs of nuclear activity (e.g., 
iRodighiero et all l2007t Ivan der Wei et all |2007) such 
that it is unlikely that central point sources affect our 
size measurements. This is also clear from the fact that 
the residuals of our R 1 / 4 profile fits generally do not show 
central point sources and that none of the deep spectra 
used to measure dispersions show evidence for nuclear 
activity. Furthermore, the good correspondence between 
the rest-frame wavelength of the imaging data sets used 
at different redshifts assures us that morphological K- 
corrections do not play a significant role. 

Despite the broad consistency between our results and 
those previously published, the agreement is not per- 
fect. There is a marginal inconsistency at the 1.5 a level 
between the size-evolution measurement from kinematic 
data and the size-evolution measurement from photomet- 
ric data shown in Fig. [6j This could point to the presence 
of some systematic effects within the z > 1.5 results. Al- 
ternatively, the different studies sample galaxies with a 
wide range in masses, and therefore mass-dependent size 
evolution could lead to apparent discrepancies among the 
samples. This is explored in the following section. 

Our robust results strengthen the results from pre- 
vious studies. We conclude that early-type galaxies at 
z = 1 are ~ 2 times smaller than local early types with 
the same mass, and that at z — 2 — 2.5 this size dif- 
ference is likely i ncreased to a fact o r of ~4, a s prev i- 
ously observed bv IZirm et all (|2007l ). iToft et alJ (|2007f ). 
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Ivan Dokkum et all (|2008D . and iBuitrago et all (|2008f) . 

6.2. Comparison with Model Predictions 

The fact that we see considerable evolution in galaxy 
size with redshift is not surprising from a theoretical per- 
spective. Most semianalytic models of galaxy formation 
in a ACDM universe predict substantial size evolution 
over the past several billion years. A comparison be- 
tween the observed and model-predicted amount of size 
evolution will help to identify the mechanism(s) that are 
responsible. In Fig. Owe directly compare the observed 
evolution in surface dens ity with the pre d ictions from 
the semianalytic work by iKhochfar fe Silk! (|2006af ) . For 
galaxies with a given mass the model significantly undcr- 
predicts the evolution in size and density, except, per- 
haps, for the most massive galaxies. In our data set we 
see no indication that the magnitude of size and den- 
sity evolution increases with galaxy mass, as predicted 
by the models. In fact, the most massive galaxies in our 
sample are precisely the only ones that are not different 
from local massive galaxies. Note, however, that statisti- 
cally speaking the evidence for mass-dependent evolution 
is weak (see ? 15. 2p . Moreover, the most massive galax- 
ies in our distant sample are a special subset, brightest 
cluster galaxies. Such galaxies have been shown to have 
properties th at deviate from those of other massive galax- 
ies (s ee, e.g.. Ivon der Linden et alj 120071 : iBernardi et al.1 
I2007T) . 

By including the z — 1.5 — 2.5 photometric samples 
discussed in § !6.1l we can place further constraints on the 
models. In Fig. [7] we compare the observed size evolution 
of the available sample s, normalized to z = 1, with the 
model predictions fromlKhochfa r~fe Silkl (|2006af) . Repre- 
senting the model predictions by a single line is justified 
by the fact that the predicted evolution of log(i? c ff) with 
log(l + z) is very close to linear. Again, the observed size 
evolution is stronger than that predicted by the model. 

It is interesting to note that, in qualitative agreement 
with the model prediction, we see a hint that size evo- 
lution depends on mass in the compilation presented in 
Fig. [7] The samples containing on average the lowest- 
mass galaxies display marginally less evolution. It has 
to be kept in mind, however, that small sample sizes 
and systematic effects are more important for determin- 
ing second-order effects such as mass dependence (see 
§ 15. 21) . A clue that systematic uncertainties may play 
a role is the remaining difference between the kinematic 
and photometric samples. Alternatively, it may signify 
non-linear evolution of log(i? c ff) with log(l + z). 

6.3. Size Evolution of Individual Galaxies 

It appears that the observed size evolution of a fac- 
tor of ^2 between z = 1 and the present for early-type 
galaxies with masses ~ 10 11 Mq is similar to the pre- 
dicted evolution for early-type galaxies that are an order 
of magnitude more massive (see Figs. [5] and [7]). This 
suggests that the mechanism responsible for increasing 
the average size of early-type galaxies with time may be 
well understood, but that it is not imple mented correctly 
in the c urrent semianalytic model from IKhochfar fc Silkl 
lj2006aH . The process of size evolution may occur at dif- 
ferent times and under different circumstances than is 
now assumed. This may be related to the late assem- 
bly of very massive galaxies in models of this kind (see 
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Fig. 6. — Size evolution with redshift as derived in this paper 
with dynamically determined masses {large filled circles) compared 
with previous results based on photometric masses (small filled 
circles). The solid line connects our samples at z ~ 0.06 and 
Z ~ 1, the dashed line is a linear least-squares fit to the small 
filled data points. The open circles are samples of cluster galaxies 
with photometrically measured masses and serve as an illustration 
that size evolution shows a continuous trend between z = 2.5 and 
the present. The broad agreement in size evolution as derived 
from galaxies with dynamically and photometrically determines 
masses reinforces the conclusions of previous, photometric studies 
whose results were potentially mitigated by considerable systematic 
effects that do not affect our analysis. 



also, e.g. jDe Lucia et al"1l2006h . a pr ediction that is chal- 
lenged by various observations (e.g.jCimatti et al.1 120061 : 
IScarlata et aHl2007l : ICool et al.ll200a^ . 

It is beyond the scope of this paper to fully discuss 
these possible discrepancies. Instead we will explore the 
question whether the proposed physical processes respon- 
sible for size evolution are consistent with the observed 
trends. In the semi-analytic models it is assumed (and 
this is confirmed by numerical simulations) that mergers 
drive size evolution. The gas content of merging galaxies 
largely determines the relative size of the merger rem- 
nant compared to its ancestors. Because gas fractions 
were higher in the past, galaxies that form early will be 
smaller than galaxies that form late. In the framework 
of cosmological simulations this means that galaxies at 
high redshift will be smaller because they were formed 
through gas-rich mergers and that those merger rem- 
nants can grow over time through subsequent mergers 
with other galaxies that are progressively more devoid of 
gas. 

The question is whether the observed size evolu- 
tion is dominated by size evolution of individual galax- 
ies or simply by the addition of larger galaxies over 
time. At z = 1 only about 30-50% of th e present- 
day early-type galaxy evolution had formed dBell et al 



2004 iBrownet aD20b7l : IScarlata et"aTll2007t iFaber et al 



2007). If we assume that these galaxies will make up the 



30-50% most dense early-type galaxies in the present-day 
universe, then the scatter in the local Md yn -^?eff relation 
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Fig. 7. — Size evolution per unit redshift vs. mean galaxy mass 
of our sample (large circle) and samples taken from the litera- 
ture (small squares; see Fig. [5] for references). Samples consisting 
of high-mass galaxies show somewhat stronger size evolution than 
samples consisting of low-mass galaxies, which is qualitatively, but 
not quantitatively, c onsistent with the predict ions from the semi- 
analytic model from Khochfar & Silk (2006a) (solid line). This 
conclusion should be considered highly tentative, however, as this 
interpretation is hampered by systematic uncertainties and small 
sample sizes. 



implies that the sample-averaged size increases by a fac- 
tor of 1.3—1.4 between z = 1 and the present. Such evo- 
lution is thus expected in the absence of size evolution of 
individual galaxies and this is less than the observed evo- 
lution of a factor of 2. To explain the observed evolution 
by growth of the early-type galaxy population without 
changes in the sizes of individual galaxies, the number 
density of early-type galaxies is required to increase by 
an order of magnitude between z = 1 and the present. 
Such strong evolution is clearly ruled out by the above- 
mentioned determinations of the number density of red 
galaxies at z ~ 1. 

Similarly, at z — 2 only ~10% of the g alaxies with 
mass es > 10 11 M had been assembled (Kriek et al. 
2008); if those galaxies, evolve into the 10% most dense 
present-day early-type galaxies then an increase in av- 
erage size by a factor of ~2 can be accounted for, less 
than the observed amount of evolution. Th ese arguments 
are in agreement with the conclusions from lCimatti et al.l 
(2008), who show that local galaxies with the same sizes 
and masses as galaxies in the z = 1 — 2 samples are so 
rare in the local universe that it can be confidently ruled 
out that their structure remains unchanged up until the 
present day. Note, however, that these arguments may 
be affected by the aforementioned biases in the SDSS 

(SHU). 

We conclude that size evolution due to the addition 
of larger galaxies over time contributes at most half of 
the observed evolution in the M^yn-RoS relation. The 
remainder must be due to size evolution of individ- 
ual galaxies. Numerical simulations have demonstrated 



that when early-type galaxies accrete neighbors without 
significant dissipational processes a c s does not change 
by much and that, to first order, i? c ff increases lin- 
early with mass. This does not depend strongly on 
the mass o f the accreted object, i.e., the mass ratio of 
the merger (iBovlan-Kolchin et al . 2005; Robe rtson et al.l 
[2^M[Bmdan-Kolchin et al.ll2006l) . 

Simulations in a cosmological context show that an 
increase in size by a factor of 2 between z ~ 1 
and the present is certainly possible (|Naab et al.ll2007l) . 
The strong observed size evolution thus argues in fa- 
vor of a scenario in which significant mass from low- 
mass companions is accreted onto existing early-type 
galaxies over the past ~7 Gyr, which also explains the 
broad tidal features that are frequently observed around 
early-type galaxies (|van Dokkuml l2005h . As shown by 
iFeldmann et al.l (|2008l ) such features are not necessarily, 
and are even quite unlikely to be, the result of major 
merger events and are most likely due to the accretion of 
low-mass, gas-poor satellites. 

We note that the size evolution of individual galaxies 
and the evolution of the sample average are insepara- 
ble because galaxies evolve in mass as well as in size. 
Nonetheless, it is important to distinguish this complex 
scenario from the simple picture in which early-type 
galaxies that form at different redshifts have different 
sizes but do not structurally evolve at later times. The 
strong observed size evolution clearly rules out the lat- 
ter, indicating that the build-up of the early-type galaxy 
population is a complex and ongoing process. 

Finally, it is remarkable that the change in the sizes 
of early-type galaxies is consistent with and differs by 
less than 15% from the change in the scale factor of the 
universe, 1 + z. Within the standard cold dark matter 
scenario this is likely a coincidence since dissipational, 
strongly non-linear processes that are decoupled from 
cosmic expansion dominate at the kiloparsec scale of 
forming galaxies. Nonetheless, we cannot exclude the 
possibility that there is an underlying, fundamental rea- 
son that galaxies are scale-invariant with respect to a co- 
moving coordinate system. In an alternative description 
of dark matter, i.e., Bose-Einstein condensed, ultra-light 
particles with a ~10 kpc-s i zed wave funct ion (fuzzy dark 
matter or FDM. [Shl[l99l iHu et alM OQO). sizes of halos 
and their occu pying galax ies possibly follow the cosmic 
expansion rate (|Leej|2008T) . 

7. SUMMARY 

In § [2] we construct a large sample of nearby (0.04 < 
z < 0.08) early-type galaxies extracted from the SDSS 
(DR6). We use the pipeline velocity dispersion measure- 
ments and obtain our own size measurements in order to 
construct the local dynamical mass-size relation (§ 12. 2p . 
In addition, we construct a sample of 50 morphologi- 
cally selected early-type galaxies in the redshift range 
0.8 < z < 1.2 with measured velocity dispersions (§ [3]). 
Sizes are determined from ACS imaging in the same man- 
ner as for the galaxies in the nearby sample, and sys- 
tematic effects are quantified through simulations (§ HJ. 
The distant sample contains galaxies in the mass range 
3 x 10 10 M Q < M < 10 12 M Q , with a typical mass of 
2 x 10 11 M Q . 

The main result is that the a e s-R c s distributions of 
the nearby and distant samples are significantly different, 
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even after we correct for the incompleteness of the distant 
sample at low masses (§EJ). The implied size evolution is 
R cS oc (l + z)-°- 98±o n , or a factor of 1.97±0.15 between 
z = 1 and the present. Similarly, the projected surface 
densities of the distant early-type galaxies are a factor of 
~4 higher than those of their local counterparts. The 
stellar populations of early-type galaxies that already 
existed at z = 1 may, for the most part, be passively 
evolving over the past 7-8 Gyr, however, their structural 
properties undergo substantial changes over that period. 

Our results are in broad agreement (see § 16. 1| with 
previously published size-evolution measurements that 
are based on samples without dynamical mass measure- 
ments and, in some cases, without spectroscopic red- 
shifts, high- resolution HST imaging, and/or consistently 
determined sizes. We therefore conclude that system- 
atic effects, most notably those in the mass estimates, 
which potentially could have hampered previous studies 
are small relative to the observed amount of evolution. 

The observed size evolution is in qualitative agreement 
with predictions from recent semianalytic models. How- 
ever, the predicted evolution is much slower than the 
observed evolution. The observed size evolution of early- 
type galaxies can be understood within the context of 
the cold dark matter scenario in which galaxies that form 
late have larger sizes than galaxies that form early, due 
to lower gas fractions at late times, and the growth of 
individual galaxies through the mostly dissipationless ac- 
cretion of satellites at later evolutionary stages. 
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TABLE 1 

Velocity Dispersions and Sizes of the 
Distant Sample 



ID 


z 


°"cff 




Reff 


MS 1 0^4 1 fS4Q 


0.831 


243 


± 


■ 


4.91 


MS 1 0^4 940Q 


0.831 


287 


± 


33 


3.30 


MCI 1 0^4 ^n^s 


0.831 


303 


± 


33 


10.20 


MS 1 f)f54-S7fi8 


0.831 


222 


± 


24 


3.28 


ms i 054-391 n 


0.831 


295 


± 


12 


1.80 


MS 1 0^4 434^ 


0.831 


336 


± 


34 


4.35 


MS 1 054-4520 


0.831 


322 


± 


30 


15.20 


MS 1 n^/L 470^ 


0.831 


253 


± 


36 


8.84 


MS 1 f)'S4-4Q2(S 


0.831 


310 


± 


38 


2.04 


MS 1 054-5280 


0.831 


259 


± 


31 


3.68 


MS 1 fl^zL ^908 


0.831 


284 


± 


39 


3.54 


MS 1 0^4-5^47 


0.831 


254 


± 


24 


2.94 


MS 1 054-5450 


0.831 


234 


± 


26 


8.16 


MS 1 054-5529 


0.831 


182 


± 


23 


3.24 


MS 1 0^4 ^77 


0.831 


305 


± 


10 


2.67 


MS 1 0^4 ^ifififi 


0.831 


286 


± 


23 


4.99 


MS 1 f)f54-fi7f5fi 


0.831 


232 


± 


27 


3.98 


ms i 054-6036 


0.831 


254 


± 


22 


2.93 


MS 1 054-6301 


0.831 


249 


± 


24 


3.55 


ms 1 054-6688 


0.831 


274 


± 


37 


2.93 


HFjFN 90fi 

L 1 LJ J_ il aUU 


0.936 


199 


± 


18 


1.11 


HDFN-237 


0.851 


280 


± 


21 


1.80 


HDFN-256 


0.974 


306 


± 


14 


3.06 


HDFN-635 


0.820 


201 


± 


17 


2.29 


HDFN 681 


0.842 


341 


± 


30 


1.43 




1.013 


374 


± 


39 


3.77 


HDFN SI 1 

i. ± _L/ J. 1 N O J- _L 


0.848 


216 


± 


14 


1.32 


HDFN-933 


0.847 


305 


± 


37 


1.86 




0.854 


235 


± 


17 


2.15 


HDFN-1 236 


0.850 


217 


± 


12 


1.98 


i. ± _L/ J_ li liiOU 


0.846 


247 


± 


17 


3.66 


L 1 _L/ J_ 1 N J. £i O I 


0.846 


342 


± 


23 


3.94 


NDFN 1 ^98 


0.845 


250 


± 


31 


1.82 


NDFN 1 ^4^ 


0.849 


280 


± 


13 


1.52 


HFjFN 1 

L 1 LJ J_ iN -LtJUc/ 


0.943 


178 


± 


12 


1.66 


L 1 _L/ J. IN 1UJO 


0.841 


330 


± 


13 


1.96 


HDFN-1 706 


0.913 


215 


± 


12 


2.23 


HDFN-1 7f)Q 

1_±-L/J_ li"l 1 UJ 


0.842 


218 


± 


11 


1.16 


CDFS-1 


1.089 


231 


± 


16 


2.83 


CDFS-2 


0.964 


200 


± 


10 


2.30 


CDFS-3 


1.044 


300 


± 


32 


1.00 


CDFS-4 


0.964 


336 


± 


19 


6.84 


CDFS-7 


1.135 


232 


± 


20 


5.77 


CDFS-12 


1.123 


262 


± 


21 


0.94 


CDFS-13 


0.980 


247 


± 


11 


2.20 


CDFS-14 


0.984 


197 


± 


23 


2.80 


CDFS-18 


1.096 


324 


± 


36 


3.97 


CDFS-20 


1.022 


199 


± 


16 


2.24 


CDFS-25 


0.967 


258 


± 


19 


0.86 


CDFS-29 


1.128 


221 


± 


18 


1.59 



Note. — The I Ds and veloc i ty dis - 
pcrsions arc taken from iWuvts et "all fe004ft, 
iTreu et all d2005bfl . and Ivan der Wei et all 
J2005I) for MS 1054-0321 , HDFN, and CDFS, 
respectively. The dispersions are all aperture- 
corrected according to eq. The sizes are 
determined as described in § 13.11 and a stan- 
dard error of 14% is adopted for all galaxies, 
which is based on the simulations described 
in §H 



